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Abstract 

We develop a systematic procedure for constructing quantum many-body prob- 
lems whose spectrum can be partially or totally computed by purely algebraic 
means. The exactly-solvable models include rational and hyperbolic potentials 
related to root systems, in some cases with an additional external field. The quasi- 
exactly solvable models can be considered as deformations of the previous ones 
which share their algebraic character. 
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1 Introduction 



Since the pioneering works of Calogero and Sutherland, [p], §, |3|, much effort has been 
invested in the study of quantum many-body problems. A thorough classification of 
completely integrable Hamiltonians related to root systems has been performed in the 
early eighties by Olshanetsky and Perelomov, both in the classical, @], and the quantum 
cases, j5|. The complete integrability of these models is associated to an underlying 
root system structure, the integrals of motion being related to the radial parts of the 
Laplace-Beltrami operator on a symmetric space with the given root system. 

Some years ago, Turbiner showed using Lie-algebraic techniques that many (though 
not all) of Olshanetsky-Perelomov's (OP) Hamiltonians are also exactly solvable, || 
By exact solvability we mean here that the iV-body Hamiltonian preserves an infinite 
increasing sequence of subspaces of known smooth functions, whereas a quasi-exactly 
solvable Hamiltonian just preserves a single finite-dimensional subspace. The construction 
of these invariant subspaces is usually based on the theory of representations of Lie 
algebras of differential operators, ||. 

The idea of constructing exactly solvable many-body problems based on the zeros 
and poles of special solutions of partial differential equations has long been known. The 
original work |J explores in depth this relation at the classical level. This idea has 
been exploited by Hou and Shifman in a recent paper, ||10|| , which extends this approach 
to the quantum case using Lie-algebraic techniques, thereby obtaining a new family of 
quasi-exactly solvable many-body potentials. 

In this paper we apply Calogero's procedure for constructing solvable many-body 
problems to the most general quasi-exactly solvable operator on the line admitting square- 
integrable eigenfunctions, finding several families of exactly and quasi-exactly solvable 
many-body problems on the line. The exactly solvable Hamiltonians include Calogero's 
original model and many of Olshanetsky-Perelomov's non-periodic Hamiltonians, and 
also a family of potentials which are not directly related to a root system due to the 
presence of an external field. The latter model has been studied by Inozemtsev and 
Meshcheryakov, [pj| , who derived its discrete spectrum from that of the standard hyper- 
bolic BCn model by a limiting process. 

It should be emphasized that, although the construction presented here leads to many 
previously known potentials based on root systems, our approach is entirely different from 
the usual one, in that no use needs to be made of any underlying root system structure. 
The main advantage of this algebraic method over the root system approach is that, while 
the latter is rather rigid, the former's flexibility allows for deformations of the exactly 
solvable models preserving to some extent their integrability properties. 

The paper is organized as follows. In Section |2| some basic definitions are given, and 
the construction of the general many-body Hamiltonian is explained. The algebraization 
of this Hamiltonian is discussed in Section |3] for the five families of normalizable quasi- 
exactly solvable operators on the line. In Section ^ the energy spectrum of all exactly- 
solvable models obtained in the previous Section eis given explicitly. In Section |5] we work 
out a few concrete examples, and we sum up the conclusions and outline future work in 
Section || Some useful expressions for handling symmetric variables can be found in the 
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Appendix. 



2 Derivation of Hn 

A Schrodinger operator 

H = -J2 d l + V ^ (1) 

k 

(where x belongs to an open subset of Euclidean space) is said to be quasi- exactly solvable 
(QES) if it preserves a known finite-dimensional subspace M. of smooth functions. The 
spectral problem for H reduces in this space to diagonalizing the matrix of H\m, which 
makes it possible to compute a finite subset of the spectrum of H by purely algebraic 
means. If H preserves an infinite increasing sequence of known finite-dimensional sub- 
spaces Aio C A4i C - • • C M.k C . . . then we shall call it exactly solvable (ES), since 
in this case one can algebraically compute an arbitrary number of eigenfunctions and 
eigenvalues of H by restricting it to each Aik- (If M. or the A^^'s are not in L 2 , some of 
the eigenfunctions of H computed in this way may turn out not to be square-integrable; 
see |12| for an analysis of this issue in the one-dimensional case.) 

The QES (or ES) character of a Schrodinger operator is invariant under a natural 
group of equivalence transformations, generated by changes of coordinates xhz = Z(x) 
and conjugation by arbitrary non-negative functions /i(z), which map H = H(x) into the 
differential operator (not necessarily of Schrodinger type) 

F(z) = M^tf (x) Mz) • (2) 

Thus, if H(x) preserves A4 = A^(x) then H(z) preserves the subspace 

A4(z) = //(z)-\M(x) , (3) 

and if ^(x) is an eigenfunction of H(x) with eigenvalue E belonging to A4(x) then 
■0(z) = /i(z) -1 -?/>(x) is an eigenfunction of H(z) lying in A4(z) with the same energy 
E. Note, however, that such an equivalence transformation may not preserve the square 
integrability of the eigenfunctions, since /x(z) is not required to be unimodular. 

A very general way of constructing QES Schrodinger operators is to start with a 
finite-dimensional Lie algebra of differential operators g admitting finite-dimensional rep- 
resentations in the space of smooth functions in the variable z (usually called a quasi- 
exactly solvable algebra in the literature). Any differential operator H(z) belonging 
to the enveloping algebra of g will automatically preserve the carrier space A^(z) of 
such a finite-dimensional representation. If one can find an equivalence transformation 
(0) such that the differential operator H(x) is of Schrodinger type (0), then H(x) is 
clearly a QES Schrodinger operator, since it preserves the finite-dimensional subspace 
A'l(x) = n{z)M.{z). In this case, one says that g = /i(z) • g ■ /u(z) _1 is the hidden 
symmetry algebra responsible for the QES character of H(x). In this paper we shall be 
almost exclusively concerned with this special type of quasi-exactly solvable Schrodinger 
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operators, that we shall call algebraic to single them out from the rest. The function /i(z) 
is usually called the gauge factor in the literature, and H(z) is referred to as the gauge 
Hamiltonian. 

In one dimension, the only Lie algebra of first-order differential operators is (up to 
equivalence) the standard projective realization of sl(2) (or its subalgebras) , with basis 
elements 

J N = d z , J° N = zd z -^-, J+ = z 2 d z -Nz. (4) 

If iV is a non-negative integer, the latter algebra admits a (N + l)-dimensional represen- 
tation in the space Vn of polynomials of degree < N. The most general second-order 
differential operator belonging to the enveloping algebra of the Lie algebra (|j), obtained 
by constructing an arbitrary quadratic combination of the generators (f|), is of the form 

-H(z) = P{z) dl + Q(z) 3 Z + R{z) (5) 

with 

Q(z) = Q(z) - P'(z) , R(z) = J2 - f Q'{z) + ^(N- l)P"{z) , (6) 

where P, Q and R are arbitrary polynomials of degree 4, 2 and 0, respectively, and the 
minus sign is for later convenience. If P is positive, the operator (|5]) is equivalent (in the 
sense of @) to a Schrodinger operator H(x) = —d% + V(x) by the change of variables 

* = / (7) 

and conjugation by the gauge factor 

H(z) = P(z)" 1/4 exp 

Let N > 1, and assume, for the sake of simplicity, that the gauge Hamiltonian (g) is 
diagonalizable in Vn- Then H(z) has N + 1 algebraically computable polynomial eigen- 
functions <pk(z) (0 < k < N) of degree < N, and therefore H(x) has N + 1 algebraically 
computable eigenfunctions of the form 

ipk(x) = n(z) ip k (z)\ z= £-i( x) , 0<k<N. (9) 

Following Calogero's original idea, consider now the time-dependent Schrodinger equa- 
tion with Hamiltonian H(x), namely 

H(x)*(x,t)=idtH!(x,t). (10) 

Since H(x) is time- independent, the latter equation will admit solutions of the form 

JV 




^(x,t) = Y,CkMx)e- iEkt } (11) 



k=0 
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where i/i k (x) is an eigenfunction of H(x) of the form (|) with energy E k , and c , . . . ,cn 
are arbitrary complex constants. Equivalently, 

${x,t)= n{z)${z,t)\ z =t-i {x) , (12) 

with 

ffMj^^t^e-^. (13) 

fc=0 

Since each (/^(z) is a polynomial of degree < N (and at least one of them is of degree 
N, since otherwise H(z) would not be diagonalizable in Vn), it follows from the latter 
equation that <&(z,t) is a polynomial of degree N in z with time-dependent coefficients. 
Therefore we can write 

N 

*(z,t) = C(t)]\[z-z j (t)]- ( 14 ) 
3=1 

It follows from Eq. (0) that each zero Zj(t) of $(<z, t) yields a zero 

xj(t) = £( Zj (t)) , l<j<N, (15) 

of the time-dependent wavefunction t). We will study the motion of these zeros, 
showing that it can be derived from a classical Lagrangian. 

Indeed, substituting fll2 ) into Schrodinger's equation ([II]) and using (|14]) we easily 
obtain the equation 

f m + f 9M + R iz)=i T m -jgco ( i6) 

^ (* - % (t)) (* - **(*)) + £ * - m + K[z) jri z - Zj (t) c(t) ' [ b) 

which must hold identically in z and t. Equating the residue of both sides at z = Zk(t) 
we arrive at the following system of differential equations for the functions z k (t): 

iz k = Q(z k ) + 2j2— { ^- = F k (z), l<k<N. (17) 

. -, Zk Zj 

Conversely, it can be shown without difficulty that the latter equations imply Eq. (|16|). In 
order to show that (xi(t), . . . , x^{t)) follows a trajectory of a certain Lagrangian system, 
we shall make use of the following Lemma: 

Lemma 1. Consider the autonomous system of first- order ordinary differential equations 
ix fc = / fc (x); k = l,...,N, x = (x u . . . , x N ) G C N . (18) 
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If the one-form Y^,k=i /fc( x ) ^ x k ^ s closed, then the second-order system obtained by dif- 
ferentiating fll8|) once with respect to t is Lagrangian, with Lagrangian given by 

N 

L = ± 2 -Y,f»- ( 19 ) 

3=1 

Proof. The second-order system obtained by differentiating (|18| ) with respect to t is 

^-iE^—EA^ l<k<N. (20) 

fe=l ■? fc=l 3 

On the other hand, the Euler-Lagrange equations associated to the Lagrangian (O) are 



N 



2 x k = — = -2 2^ /j 



— ^, f < fc < N. 

dx k 



(21) 



If X^fcLi /fc( x ) is closed then 



df k dfj 



dxj dxk 



so that (12(1) and (ETJ) are indeed identical. 



Q.E.D. 

In our case, the system (|18|) is obtained from ( |17D by the change of independent variables 
(|T5|), with ^(z) as in (0). Since Xk = P(z k )~ 1 ^ 2 z k} we have 



f k = p{z k y l/2 F k 



and 



N 



N 



Fdz) 



fc=i fc=i v K > 

The latter one-form is clearly closed, since if j ^ k we have 



(22) 



(23) 



d fF k (z) 



dzj \P(z k ) J dzj 







Q{zk) 
P(z k ) 



A' 



+£ 

z=i 
i^k 



Zk — z i 



[Z k - Zj) 2 



which is clearly symmetric under the exchange of j and k. It follows from the previous 
Lemma and Eqs. ([19|) and ( p2|) that the zeros x(t) of the time- dependent wavefunction 



11|) move along a trajectory of the Lagrangian system with Lagrangian 



t _fv f 



(24) 
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whose associated Hamiltonian is 

N N T? 2 ( \ 

k=i k=i ^ k > 

Substituting ([H]) into ( j25| ) and letting p k = —id Xh we obtain the iV-particle quantum 
Hamiltonian 

N 

H N = -J2dl+V N (x), (26) 



where 



k=l 



N A 2 N ^ N „ AT 



fc=i rk j,k=i Zk i j,k=i Z kj j, k ,i=i Zk i Zkl 



and we have set 

P k = P(z k ), Qk=Q(z k ), Z k j = Z k -Zj. 

It will be convenient for our purposes to allow the coefficients of each of the four sums 
in Eq. (|27|) to be different. Dropping an inessential overall factor we finally arrive at the 
following formula for the potential: 

i N r> 2 N n N p N p 

*w<> = iE%+»E£+*Ej+*E (28) 

k=l ft j,k=l K3 j,k=l k 3 j,k,l=i 3 

The main purpose of this paper is to show that, for certain choices of the polynomials 
P(z) and Q(z) (and, sometimes, the constant gi), the iV-particle quantum Hamiltonian 
with potential (|28|) is quasi-exactly solvable or even, in certain cases, exactly solvable. 
This will be proved in the next Section, essentially by showing that Hn lies in the en- 
veloping algebra of a quasi-exactly solvable Lie algebra of first-order differential operators 
equivalent to a certain standard realization of sl(N + 1) (cf. Eq. fl29|) below). 



3 Algebraization of Hjy 

We shall prove in this Section that the Hamiltonian (j26"D-(f28|) derived in the previous 
Section is algebraic, provided that P(z), Q(z) and (in certain cases) g\ are chosen appro- 
priately. In analogy with the one-dimensional case, we shall show that H^(pc) is equivalent 
under a change of variables and a gauge transformation fl2|) to a gauge Hamiltonian Hn(z) 
which can be written as a quadratic polynomial in the differential operators 

V k = d Tk , Mjh = Tj d Tk , U k = r k - E n d n 
j,k = l,2,...,N, (29) 
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spanning the Lie algebra sl(N + 1). In the previous formulas m is a non-negative integer, 
and we have denoted by 



h<h<— <ik 



the /c-th elementary symmetric function. Since the operators ( p9|) preserve the finite- 
dimensional module 



N 



M m = span jrf 1 ^ 2 • • •rjf : J^h < m\ (31) 

of polynomials of degree < m in the symmetric variables r = (n, . . . , Tjv), H^(z) is 
algebraic, possessing (at most) dim.M m = ( m ^ 7V ) eigenfunctions Xfe('7"( z )) that are poly- 
nomials of degree < m in r. It follows from the discussion in the previous Section that 
ifjv(x) is also algebraic, and that it admits eigenfunctions of the form 

V fc (x) = ^(z)xfc(r(z)), (32) 

where /z(z) is the gauge factor, whose specific form will be given below (cf. (|37|)), and 
X/c(t(z)) is a polynomial eigenfunction of ifjv(z) which can be algebraically computed. 
For this reason, we shall henceforth use the term algebraic eigenfunction to refer to an 
eigenfunction of of the form (j32|). 
The generators 

T> k , A/}*; l<j,k<N (33) 

span the Borel subalgebra b^+i C sl(iV + 1), which preserves the infinite sequence of 
subspaces .Mo C M.% C .M2 C . . . . Therefore, if Hn is a quadratic combination of the 
differential operators ([33]) then the physical Hamiltonian Hn is exactly solvable. 
It will also prove convenient in what follows to define the polynomial subspaces 

C k = span jr^ 2 . . . rjf : ^ 2 = k\ (34) 

of all the symmetric polynomials homogeneous of degree k in the variables z = (z±, . . . , zn), 
and their direct sums 



M m = ($£ k . (35) 



fc=0 



Clearly, Xi m C -M m , though in general jti m need not be invariant under the action of 
Hn. An important exception occurs when is exactly solvable; indeed, it will be shown 
in Section [| that in this case Hn preserves the infinite sequence .M0C.M1C.M2C.... 
This important fact shall be used in Section [| to derive a formula for the energy spectrum 



8 



of all the exactly solvable models we shall construct, and in Section [5] to exactly compute 
some eigenfunctions of and their corresponding energies for an arbitrary number of 
particles N. 

We shall restrict ourselves in this paper to polynomials P(z) whose corresponding one- 
dimensional QES operators (H) admit normalizable (i.e, square-integrable) eigenfunctions, 



12], leaving the periodic cases for a future work. In this case, there are 5 canonical forms 



for P(z), all of which can be taken as monic polynomials of degree not greater than 2. 
By (|6p, Q(z) is then an arbitrary quadratic polynomial 

Q(z) = c + z 2 + c z + 5- . (36) 

Let us define the gauge Hamiltonian H N (z) by Eq. (Q), where the gauge factor 



Mz)=n4-n p i ex p(| r 

j<k k I 



(37) 



has been chosen by analogy with the one-dimensional case (cf. fL|), and a, 6, c are real 
parameters. (From now on, all indices in sums and products will implicitly run from 1 to 
N, with the restrictions indicated under the summation or product symbols). Using (p8|) 
and dropping an additive constant we obtain the following explicit formula for ifjv(z): 



H N (z) 
where 



u u L \ / 



&+2 ) P 'k + cQk 



V N (z) = A 1 r 1 + A 2 J2^ + J2y 

j+k k 3 k k 

and the coefficients Ai are given by 



A 3 pf + A 4 Q k P^ + A 5 Ql 



(3J 



(39) 



A x = c + [{N-l)( gi -ac) -c) 
c 



A 4 



4 



(2b 



A 2 
A 5 



g 2 - a(a - V 



4 



.4, 



(40) 



We want to find sufficient conditions for to be expressible as a quadratic combi- 
nation of the generators (p9|). From (^) it is clear that one such condition is A 2 = 0, 
which yields the following relation between the coupling constant g 2 and the exponent a 
in the gauge factor (|37j): 



g 2 = a (a - 1] 



The differential part of if at equals 



— m cc + Ti + c c + Ui mod b jv+i 

= —m c c + 7~i mod sl(N + 1) . 



(41) 



(42) 
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In Cases 1-3, P(z) = z 2 + e with e = 0, ±1, and the last term in (|39|) can be expressed 
(up to a constant) as 

A, cl + 2~c + (A, + ~c A 5 ) r l + J2 ■ (43) 



k 



for certain constants po an d P\- It is then clear that the terms proportional to T\ in (p9[) 
and (H2|) must cancel, and the remaining terms in (S) must vanish, giving rise to the 
conditions 

A 5 c+ = (44) 
Po = Pi = (45) 



(-Y- !)((/, -(/r) - (6 + m +i) r 



. (46) 



The analogous conditions for Cases 4 and 5 (when P(z) = z and P(z) = 1, respectively) 
can be found at the end of this Section, when we discuss those cases in detail. 

In particular, if c + = then the "quantization condition" fl46|) is automatically sat- 
isfied and the term proportional to the generator U\ in (|42|) vanishes (this also holds in 
Cases 4 and 5), so that in this case is a quadratic combination of operators belonging 
to the Borel subalgebra of sl(N +1). Therefore the condition c + = ensures the exact 
solvability of the potential (p8|) . The specific conditions for H at to be algebraic, together 
with the explicit form of the potential ( p8|) and the normalizability conditions for the 
algebraic eigenf unctions, will be discussed on a case by case basis below. 

All the potentials discussed in what follows diverge when Xj = Xk (cf. Eq. (^)). For 
this reason, we shall consider them to be defined in the open set 

xn < xn-i < ■ ■ ■ < xi . (47) 

Moreover, the finiteness of the mean kinetic energy of the algebraic eigenfunctions near 
the hyperplanes Xj = Xk requires that a > 1/2, so that g 2 > — \ in all cases. In Cases 2 
and 4 we shall impose the additional restriction Xk > (1 < k < N) to make the change 
of variables (0) one-to-one ((|28|) is generically singular at Xk = in these cases), so that 
the potential is defined in the open subset 

< x N < x N _i < ■ ■ ■ < x x . (48) 

This corresponds to choosing a fundamental chamber for the action of the Weyl group, 
I- 

The following constants will be used throughout this Section: 

C± = cc±, C = c c , (49) 

B = b+-cco, B_ = b + cc_. (50) 
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3.1 Case 1. = z\ + 1 , = sinhx^. 



It can be shown that a necessary and sufficient condition for all the wavefunctions in the 
algebraic sector to be normalizable is 



c + = 0, 

a{N - 1) + B + m < 0; 



(51) 
(52) 



hence (pH| ) and (fR]) are automatically satisfied, and the potential in this case is exactly 
solvable. Conditions (PEBT) reduce to 



C-[B--)= m , 



(53) 
(54) 



where rj and 771 are defined by 



1 



(55) 



The solution to the previous system becomes unique after imposing the normalizability 
condition fl5"2"|), and is given by 



5 = i(l-V5) 



(56) 



with 



s _ 1 
2 ~ 4 

The potential (|28|) is in this case 



Vo + m ■ 



(57) 



Vjv(x) = '^2(Vo + Vi sinn x k) sech 2 x fc + g 2 tt 



cosh 2 Xfc 



(sinh Xj — sinh re*. 



,2 • 



(5? 



which, after dropping the constant term g 2 N(N — l)/2, can be expressed in the more 
familiar form 



V N (x) = ^2 (Vo + Vi sinhx fc ) sech 2 



Xk 



csch' 



2 / Xj Xk 



sech z 



2 f Xj ~\~ Xk 



(59) 
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In the latter formula, r] and % can be taken as arbitrary real parameters (cf. (|55|)) 
restricted only by the normalizability condition (|52"D, which can be written as 



m< - (Vs-l) -o(JV-l) 



(60) 



with s defined by Eq. (ply} . The exactly solvable potential ( p9|) will admit square- 
integrable algebraic eigenfunctions if the latter inequality can be satisfied by a non- 
negative integer m, i.e, if the right-hand side is positive, which yields the restriction 

l2 



s > [l + 2a(iV-l)] 
The algebraic eigenfunctions ( |32|) in this case have the form 

\ a 

Xj ~T~ Xk 



(61) 



^(x) = x(r) • JJ ( cosh' 
j<fc 



sinh — 



] J sech 



B| — s— arctan(sinh Xfe) 

Xfc e 



(62) 



where a, 5, and C_ are given by Eqs. (|4T|), (|56|), and (|57|). The function x(r) is a 
polynomial eigenfunction of of degree < m in the symmetric variables r^, where m 
is the largest non-negative integer satisfying (|6Tf ). 

The potential (59) is not, as it stands, one of Olshanetsky-Perelomov's completely 
integrable models. If we perform the imaginary translation Xk *> + i f (1 < k < AT), 
then fl59D becomes 



ft- (ft 



/3-(ft 



sinh 2 ^ 



cosh 2 ^ 



+ 



92 



csch 



+ csch^ 



with 



P± = B T i 



(63) 



(64) 



which is a completely integrable potential of BCn type. Note, however, that the coeffi- 
cients (3±((3± — 1) are, generally speaking, complex for real values of B and C_. In fact, 
these coefficients can be real only in two cases: i) C_ = 0, and ii) B = 1/2. In the first 
case, both parameters ft. = B are real and equal, and ( |63[ ) reduces to the well-known 
completely integrable potential of Cn type 



V N (x) = B(B ^ch 2 x k + J2 

k j^k 



csch 



Xj + Xk \ ,2 

1 + csch 



(65) 



Indeed, in this case the algebraic eigenfunctions (|62| ) reduce to 



^(x) 



X{T, 



■ Y\ ^sinh 



j<k 



— sinh 



(66) 
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thus recovering the results in Ref. ||. The second case is not interesting, since B = 1/2 
implies 

P±(£t - 1) = -\ (i + ci) . (67) 

Therefore in this case the Hamiltonian is not self-adjoint, since the potential behaves as 

-\^ + Cl)xf<-\xf (68) 



as Xk — > 0. Thus when C_ 7^ 0, i.e, when r\\ 7^ 0, the exactly solvable potential ([59]) does 
not seem to be directly related to a root system. 

3.2 Case 2. = z\ — 1 , = cosha^. 

In this case, the normalizability of the algebraic eigenf unctions requires that either 

C + < (69) 

or 

C+ = and a(N - 1) + B + m < , (70) 

which leads to two different subcases, the first one consisting of QES and the second one 
of ES potentials. We shall therefore discuss separately the two cases C + < and C + = 0. 

3.2.1 Case 2.a. C + < 0. 

Conditions (HP and (f45D reduce in this case to 



c 2 = 1 (71) 

C + + C_ = (72) 

B(fl-l) = S (73) 
while condition (^) can be expressed as 

c(N- l)^i = o(iV- 1) +m + b + i (74) 

Imposing regularity of the wavefunctions on the hyperplanes = forces a unique 
solution to Eq. fl73f) , given by 



5 = 1(1+^/1 + ^). (7r,) 
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Taking into account Eqs. ([7l| ) -([74] ) , the potential 
'CI 



V N x 



E 



4 



sinlr x k + C + a(N - 1) + B + m 



can be written as 
1 



C 2 

cosh Xk + csch 2 Xk 



csch 



2 / *£fe 



+ csch z 



(76) 



This potential is quasi- exactly solvable, and it has dim.M m = ( m ^ 7V ) algebraic eigenfunc- 
tions of the form 



V>(x) = x(r) ■ JJ f sinh Xj + — sinh Xj ^ Xk \ ■ JJ(sinhx fc ) B e 



2i C osh Xk 



(77) 



where .B is given by Eq. (|75|) , and x( r ) is a polynomial eigenfunction of H ^ of degree < m 
in the symmetric variables r^. The potential (|76|) is regular on the hyperplanes Xk = if 



and only if C$ = 0. In this case 5 = 1, and therefore the algebraic eigenfunctions (77) 
are naturally extended as odd functions to the region 

|xjv| < [rrjv— 1| < • • • < . (78) 

A potential somewhat more general than ([F6|) has been proved to be completely in- 
tegrable in the classical case by Inozemtsev and Meshcheryakov, ||13|| , fi~4}| . Its quantum 
integrability has been conjectured by the same authors in Ref. [[11], although, to the best 
of our knowledge, no proof of this fact has been given. We have shown in this Section 
that the potential (|76|), which belongs to the general class considered in Ref. [11], is quasi- 



exactly solvable. In fact, we shall explicitly compute in Section [|some eigenfunctions and 
eigenvalues of the Hamiltonian with potential (|76|) for m = 1 and N = 3 particles. 



3.2.2 Case 2.b. C+ = 0. 

In this case the potential (|28|) , which is exactly solvable, can be written as 
Vat(x1 



^ = 4 E K csch2 (if) ~ V- sech - 2 (y 



where 



Conditions ( f44[) and 



V± = \(c--±5 ) 2 >0. 



( sch 2 ( 1 + csch 2 ( Xj Xh 



are now automatically satisfied, while ([|5|) reduces to 



- (r/+ - ?7_) 



(79) 
(80) 

(81) 
(82) 
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Introducing the parameters 



the latter system reduces to 



P± = B±^, (83) 



P±(P±-l)=V±, (84) 



which uniquely determines (3±. Indeed, the algebraic eigenfunctions in this case are of 
the form 



^(x) = x(r) • Yl ( sinh %3 + ^ Xk sinh X] Xk j • JJ sinlA ) coslr 5 - ) 

j<k \ ~ / k 



(85) 



where x( r ) is a polynomial eigenfunction of H jy of degree < m in the symmetric variables 
r. Therefore (3 + > 1/2 for the eigenfunctions to be regular (with finite kinetic energy) on 
the hyperplanes Xk = 0, while /3_ < —1/2 is necessary for the normalizability condition 
(|70|), which now reads 

^ (/?+ + /?_) + a(iV-l) + m<0, (86) 
to hold. Thus Eqs. fl8HD have the unique solution 

P ± =±(l± , (87) 

from which it actually follows that (3 + > 1 and (3_ < —1. 

The exactly solvable potential ( |7S|) is known to be integrable, and is associated 
with the root system BCn, when i] + ^ r]-, or Cn, when 77+ = r/_. The potential is 
regular on the hyperplanes Xk = if and only if (3 + = 1, in which case the algebraic 
eigenfunctions (p5|) are naturally extended as odd functions to the region ([78]). 



3.3 Case 3. P k = z\ , z fc = e 3 *. 

For the algebraic eigenfunctions to be square-integrable, one of the following three con- 
ditions must be satisfied: 

1. C+ < 0, c_ > 

2. C + = , C_ > , a(iV - 1) + B + m < 

3. C+ < , C_ = , B > 

In principle, the first case yields a family of QES potentials and the remaining two cases 
a family of ES ones. However, it can be shown that the potential in the third case can 
be obtained from that in the second one by the reflection x 1— > — x. We shall therefore 
consider only the first two cases. 
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3.3.1 Case 3a. C + < , C_ > 0. 



Conditions ( f44| ) and (45) are satisfied if 



and the quantization condition (|4q) can be written as 



c </i (JV - 1) = a(iV - 1) + 6 + m + - 



The potential (p8[) becomes in this case 

Y C 2_ 



(88) 



*9) 



^(x)=^ 



^ + 1 Co C _ e -^ + C+ ( a(A r _ i) + ^2 + m + i) e ^ + 



C3 



#2 
4 



csch^ 



(90) 



The quantum Hamiltonian (|26|) with potential (|90|) defined in the region ([17]) is quasi- 
exactly solvable. It possesses ( m ^ 7V ) algebraic eigenfunctions of the form 

^(x) = x (r) ■ J] ~ eX T • n e " (Co+1K " e^ c+eXk -°- e - Xk ), (91) 

j<k k 

which can be more conveniently written as 

V,(x) = X (r) n sinha (^y^) . JJes^^+i]'* e h(c + e* k -c.e-* k )^ (92) 

where x( r ) is a polynomial eigenfunction of if at of degree < m in the symmetric variables 

The potential ( p0|) with arbitrary constants multiplying each exponential has been 
shown to be classically completely integrable in Ref . |H| , but the explicit solutions of the 
equations of motion could only be calculated when certain restrictions on the constants 
were imposed, and only for a subset of initial conditions, [|15|]. In the quantum case, 
no exact formulas for the eigenfunctions and their corresponding energies are available 
when both the negative and positive exponentials are present, although these models 
have been conjectured to be completely integrable also in this case, lfiTJ| . We have proved 
that for certain values of the coefficients multiplying the exponentials the potential (p0|) 
is quasi- exactly solvable, and thus some of its energies and eigenfunctions (in general 
corresponding to excited states) can be exactly computed. See Section |5] for a concrete 
example when m = 1. 



16 



3.3.2 Case 3b. C , = , C_ > 0. 



The conditions ( f44|) and (^) lead again to (|88|) , but now the quantization condition (f46|) 
is absent. The potential can be simply obtained by making C + = in ( fjoj) : 

^(x) = Efx e " 2 " + ^° c - e H + f D csch2 • (93) 

fc L -I j^k ^ ' 

Note that this potential depends effectively only on one parameter, since the ratio |C_/C | 
can be assigned any prescribed positive value by performing a suitable translation of the 
coordinates. The same argument shows that the discrete spectrum of the potential (|93D 
cannot depend on C_. 

The quantum Hamiltonian ([26|) with potential ( |93D is exactly solvable. It describes an 
external Morse potential acting on each particle, plus an A N -type hyperbolic interaction. 
The discrete spectrum of this potential has been computed by Inozemtsev and Meshch- 
eryakov, jLl]], by relating (|9~3D to the integrable BCn potential (|79"1) through a certain 
formal limit. 

The square-integrability of the algebraic eigenfunctions is ensured if Eq. (|5^) holds 
or, taking into account (f 



Co + 2a(iV - 1) + 2m + 1 < . (94) 

This inequality can be satisfied for non-negative integer values of m provided that 

C + 2a(iV- 1) + 1 < 0. (95) 

In particular, note that the latter inequality implies that Cq is negative, and therefore 
the coefficient of e~ Xk in the potential is negative as well. The algebraic eigenfunctions 
are given in this case by 

= X(r) ]Jsmh a (^L_^\ JJexp (hc Q + a(N - 1) + l]x k - ^e~A , 

j<k \ / b V / 



(96) 



where x{ T ) is a polynomial eigenfunction of of degree < m in r, m being the largest 
non- negative integer compatible with (|9~4|). 

The limit C_ = corresponds to one of Olshanetsky-Perelomov's integrable potentials 
associated to the root system, namely 

V N {n) = -a(a - 1) Y, ^h 2 (^p) , (97) 

whose algebraic eigenfunctions are no longer normalizable (indeed, this potential is known 
to have no bound states, |J). 
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3.4 Case 4. P k = z k , z k = 

The algebraic eigenfunctions will be square integrable provided that 

C + < 



(9* 



or 



C 4 







Cn < 0. 



(99) 



Conditions (fl4|)-(|46|) must be modified in this case (since P is no longer quadratic), but 
after applying similar considerations it can be shown that the operator H ^ is algebraic 
provided that 



c 2 = l 



c 4 



b 3 

(a-c gi ){N-l)+m+- + - 



0. 



The potential is given by 
Vn(*) = £ 

k 

where 



x k. 



Vo 

V2 
V3 



16' Vl 



C0C4 



+ 



1 



[Xj - x k ) 



(100) 
(101) 

(102) 



(103) 



and 



ci,c± 

16 4 
C 2 > 0, 



5_ 



32 ' 
a(iV - 1) +m 



(104) 



1 + 4773 . 



(105) 



(We have taken the positive square root to ensure that the algebraic eigenfunctions are 
regular when Xk — > 0.) When C + 7^ the potential ( |103| ) is quasi-exactly solvable. The 
coefficients in the potential are not all independent in this case, but satisfy the relation 



jrj2 

4 ^o 



B- 

a(N - 1) + m H h 



The algebraic eigenfunctions are of the form 



^(x)=x(r)-nK-4) a -n 



(106) 



(107) 
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where B_ is given by ( |105| ), and x( T ) is a polynomial eigenfunction of H ^ of degree < m 
in t. This QES potential has been recently obtained in somewhat less generality by Hou 
and Shifman, |T0| . 

If C + = then the potential (|103|) becomes exactly solvable, and assumes the simpler 
form 



^)=e(t^ + !) 



92 Yl 



+ 



1 



(Xj + X k ) 2 (Xj - x k ) 



(108) 



where now the parameters to = 2 v /r/2 = |Co| /2 > and r/3 > are independent. The 
algebraic eigenfunctions are in this case 



^)=x{r)-\{{x)-xir.\{ 



B- -^ x 2 

x u e * *= 



j<k 



(109) 
1, and 



where B_ and x( T ) are as before. This potential is well known to be integrable, 
is commonly known as a rational B^-type potential with harmonic force. 

The potentials ( |103| )- (|108|) are regular on the hyperplanes x^ = provided that ^3 = 0, 
or, equivalently, B- = 1. By ( |107| )-( |T0~9| ), in this case the algebraic eigenfunctions can 
again be extended as odd eigenfunctions to the region (|78D . 

Following the idea of Minzoni, Rosenbaum and Turbiner, ||16|| , we can obtain a different 
QES deformation of the ES potential ( |108| ). Note that the gauge hamiltonian in this 
case can be written as 

H N = + , 



where 



= - 
A = N 



■r l dl-{A + C Q r l )d T 



B_ 



+ a(N-l) 



(110) 
(111) 



depends only on n, and 



h® ■ 0(n) = . 

If we consider eigenfunctions x °f Hjq depending only on T\ then the problem reduces to 
the one-dimensional equation 

-r lX " -(A + ConW = E X . (112) 

If we add a function U(t\) to the potential ( |108j ), we can still look for eigenfunctions of 
H 7v depending only on n, which must satisfy the equation 



n X " + (A + CorOx' + [e- U(n)]x = o 



(113) 



A wide class of exact solutions of ( [L13| ) is obtained when the operator h = + U(ti) is 
equivalent under a gauge transformation to a quasi-exactly solvable operator on the line. 
The most general such operator equivalent to h must be of the form 



h 



ri d 2 + (2 7 t 2 + Pn + a) d T1 - 2n 7 n + E 



(114) 
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where E is a constant and n is a non-negative integer, [12|. The latter operator leaves 
invariant the space of polynomials in T\ of degree not greater than n, and consequently 
admits n + 1 algebraic eigenfunctions. For h and h to be equivalent we must have 



h = e" /(n) • h ■ e /(Tl) . 

From this relation we obtain the equations 

2n /' + A + Co n = - (2 7 r 2 + /3 r x + a) 
Ti(f" + f' 2 ) + (A + C n)f -U = 2n 1 r 1 -E , 

from which we easily get 

f = -\ [irl + (/? + C )n + (A + a) logrj 
and (after dropping a constant term) 

y(r 1 )=7 2 r 1 3 + /3 7 r 1 2 + ^/j 2 -C 2 + 4 7 (a-2n-l)]r 1 - (A ~ a ~ 2)(A + a) 

4 4 n 



If we define 



E 



(115) 

(116) 
(117) 

(118) 
(119) 



4n, 



then we have shown that the iV-particle potential 



V^x) = Vj^x) + U{r 2 /A) 



D 



Ar b + Br /l + Cr z + — + B_(B_ 



+ a(a — 1) 



+ 



(Xj- + X k ) 2 (Xj - x k ) 



with 



B 



16 



D 



-(A- a- 2)(A+a) 



possesses n + 1 algebraic eigenfunctions of the form 



^fc(x) 



-(A+a) 



7 4/3 » 

<>X1,) 1 ~32 r ~ 8 r 1 



/4) • n< 

j<k 



™2 _ 2\a 



n 



k ' 



(120) 



(121) 
(122) 

(123) 



where cr^(ri) is a polynomial eigenfunction of /i of degree less than or equal to n. These 
eigenfunctions are square-integrable and regular at the origin provided that 



7 > (or 7 = and (3 > 0), a > 



N 



B_ > 



(124) 
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3.5 Case 5. = 1 , Zk = 

The conditions for i/^r to be algebraic reduce in this case to 

A 5 c + = A 5 d = . (125) 

On the other hand, the necessary and sufficient conditions for the algebraic eigenfunctions 
to be normalizable are 

C+ = 0, C <0. (126) 

Hence the all potentials in this class are exactly solvable. From Eqs. ( |125| ) and (|126|) it 
follows that c 2 = 1, so that the potential (|39|) reduces to the celebrated Calogero model 

V N (*) = 4 + 92 £ (g ._! Xfc)2 ' ( 12? ) 

with 

IC I 

w=^>0. (128) 

(To obtain the previous formula, we have performed a constant translation to get rid of 
the irrelevant parameter C_.) The algebraic eigenfunctions are of the form 

^(x) = x (r(x)) • - **)■ ri e f xl > ( i29 ) 

where the functions x( r ) are polynomial eigenfunctions of H^, known in the literature 
as the Calogero polynomials. 

4 Energy Spectrum 

In this Section we shall obtain an explicit expression for the energy spectrum of all the ES 
models derived in the previous Section, discussing some of its basic properties. The key 
idea in this respect is to show that the ES Hamiltonians preserve not only the subspaces 
Aik, but also the smaller subspaces M.^. We shall then define a lexicographic ordering of 
the basis elements within each subspace Aik ensuring that the matrix of the restriction 
of Hn to M.k is triangular. The spectrum of H^, and therefore that of Hn, will thus 
simply consist of the diagonal elements of these matrices for k = 0,1,2, ... . 



4.1 Cases 1-3 

The gauge Hamiltonian ( [38|) corresponding to the ES Hamiltonians with potential (|59|). 



and fl93|) can be written as 



H N (z) = E 



2 



Z kj 
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where e = 1, —1, 0, respectively, and 



En 



-N 



B+ a -(N-l)]\^(N*-l)}. 



(131) 



For the potentials (|59|) and ( ff9|) 5 and C_ are given by Eqs. (|56|)-(|57D and (|83|), respec- 
tively, while £> = | (Co + 1) for the potential (|93"1) on account of fl88f) . The expression of 
Hn in terms of the symmetric variables is easily found to be 



H N = E -J2 A^d Ti d Tj - (25 + 1) ^>A - a^j(2iV - j - l) rj d Tj 



)T,-id T 



i 

+ e ( - E 4 0) ^^ + a Z> - j + - j + 2)^. ) , 



(132) 



where the coefficients can be found in the Appendix. 

It follows from the latter expression and the structure of the coefficients that H ^ 
preserves the subspaces 



Mi 



j=0 



for arbitrary k, where 



N 



Cj = span < rj 1 t 1 2 2 . . . rjlf : i U 



i=l 



In fact, the first line of ( |132| ) leaves the subspace Ck invariant, while the second and third 
lines take an element of Ck into Ck-i and £fc_2, respectively. 

We now introduce an ordering of the basis of Aik consisting of monomials. First of 
all, we declare the monomials belonging to Cj to be less than those in Ck if j < k. Within 
each Cj, the ordering is then defined as follows: if the monomial rj 1 . . . r^f is denoted by 
the multi-index I = (li, l 2 , . . . , In), then 



/ < /' if Im — I'm, I 



i+1 



''i+l' H ^ L i ■ 



(133) 



Thus, for instance, Tj < T\Tj_\ < . . . < t(~ 2 t 2 < t{. 

Given this ordering of the basis, it is straightforward to check that the matrix 
is upper triangular. Its diagonal elements 



En 



2B + a{2N - 1) + 



(134) 
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give therefore the energy spectrum of H N , and hence of H N . The spectrum consists of a 
term 

E - [2B + a(2N - 

i 

which is constant over Cj, while the remaining two terms split the energies within this 
subspace, although in general they do not remove all degeneracy. It can be proved with 
the aid of the normalizability condition that the lowest energy corresponding to each 
block Cj is Ej while the highest one is i?o...i...n where the 1 is in the p-th position, 
and j = N r + p with 1 < p < N — 1. This easily implies that the lowest energy in each 
block Cj is lower than the lowest energy in the next block Cj + \. In particular, the first 
excited state always lies in .Mi (since dim .Mo = 1) ; an d can thus be exactly computed 
for any number of particles. On the other hand, some energies in Cj + i might be lower 
than those of Cj. Note also that the formula ( |134j ) for the eigenvalues of the Hamiltonian 
with potential (|7^) and (|93|) is considerably simpler than the one given in Ref. (Tl| . 



The spectrum of the three ES potentials ([5^), (|7£|) and fl93|) is the same when expressed 
in terms of B and C_, although B and C_ are differently related to the parameters 
appearing in these potentials. Due to the triangular nature of the matrix representing 
H N in M. k) it should not be difficult to give recursive expressions for the eigenf unctions, 
though we postpone this task for a future work. 

4.2 Cases 4-5 

The spectrum of the potentials ( |108| ) and ( |127| ) is highly degenerate, since in this case the 
energy E^...^ is easily seen to depend only on the single quantum number j = J2f=i i h- 
More precisely, using the explicit expression fl38|) for the gauge Hamiltonian and passing 
to the symmetric variables it is straightforward to show that 

E 3 =E + 2coj, j = 0,l,2,..., (135) 

where the ground state energy E is given respectively by 

E^ =N[1 + 25_ + 2a{N - 1)] | , (136) 
Ef ] =N[l + a{N . (137) 

We have thus rederived (albeit making no use of the underlying root systems) the well- 
known result that the energy levels in these models are equally spaced, p[. 



5 Examples 
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5.1 ES potential from Case 1 

For the ES potential (|59|) 



V N (x) = (r/ + rji sinhx fc ) sech 2 



4 ^ 



csch 



2 / Xfc 



sech 



3j q | iC 



(138) 



the gauge Hamiltonian in the symmetric variables is explicitly given by 

j 



H 



N 



d T .d T 

1 i 1 



d T .d T . 



-E 
-E 

»>3 

-Y^j[2B + l + a(2N - j - 1)] Tj d Tj - C_ ~ J + 1)^-1^ 

+ a ^(iV - j + 1)(JV - j + 2)r,_ 2 <9 Tj + E , 



k=l 

3-1 

(N -i + ly^Tj-i - ^(z — j + 2A;)r i+fc _ir i _ A ._ 1 

fc=l 



(139) 



where 5, C_, and £?o are given by Eqs. (|5^)-([57|) and ( |131| ), respectively. We know from 
the previous Section that the gauge Hamiltonian ( |139| ) preserves each subspace Aik for 
arbitrary k. For instance, the restriction of if at to M. 2 has the following matrix with 
respect to the basis {1, r±, T2, r\ } of M. 2 '- 



where we have set 



fEo -iVC_ aN(N - 1) -2N \ 

E + l-f3 -(JV-l)CL -2NC- 

E + 2(a-p+ 1) 4 

\ E -2p/ 



^ = a(N-l)+B + l<-l. 



(140) 



(141) 



(The latter inequality is obtained by imposing the normalizability of the algebraic eigen- 
functions belonging to M. 2 .) We can exactly compute four eigenfunctions of the potential 
( |138| ) with their corresponding energies by diagonalizing the matrix ( |14U| ). The energies 
are simply the eigenvalues of this matrix, namely its diagonal elements 

Ex = E - (3 + 1 , 

E 2 = E -2(3, 

E 3 = S + 2(o-/3 + l). 
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The corresponding eigenfunctions are given by 



< i < 3 



(142) 



where 



and 



with 



// = Y\ ( cos h 

j<k 



Xo = 1 



Xj + Xk . n Xj Xfc 

— sinn — 



£[(sechx fc ) |B| e~ ^Wsinh^) 



Xi = n + 
X2 = r x 



NC. 



(3-1 

, 2r 2 2C-(Na + l) 



+ 



N{aN + l)(C 2 +(3 + 1) 



X3 = T 2 



a + 1 (a + f3(a+l)((3 + l) 

C-(l-N) N(N -l)(2a 2 -a(3 + a + C 2 _) 



2a- (3 + 1 



n + 



Ti 



sinh 



2(a-(3+ l)(2o-/3 + l) 



r 2 = sinh sinh 

j<fc 



(143) 

(144) 
(145) 

(146) 

(147) 

(148) 



Note that, since ipo = \i has no zeros in the configuration space ([|7|), it is the ground 
state of the system. We have thus been able to exactly compute some eigenvalues and 
eigenfunctions of the potential ( 138 ) for an arbitrary number of particles N. 



5.2 QES potential from Case 2 

For simplicity, we shall restrict the number of particles to N = 3 and set m = 1 to find 
4 eigenvalues of the potential ( f76"D 



%(x) = E 



k=l 



c 2 

B(B — 1) csch 2 Xk + sinh 2 Xk + C + I 2a + B + - ) cosh Xk 



a(a — 1) v-^ 

j,k=i 



< sch 2 ( 1 + csch 2 f Xj Xk 



, (149) 



which can be thought of as a deformation of a Cjv-type hyperbolic potential. The gauge 
hamiltonian H% = yT l H% // can be written as 



fe=i 



fc=i 



fc=i 
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where z k = coshx^, 



H = I j ( sinh 

i<i<fc<3 



Xj ~t~ X k . 1 Xj Xfc 

— sinn — 



J^(sinh x k ) B e~ 



■ cosh Xk 



(150) 



k=l 



and 



e = - (3B 2 + QaB + 5a 2 ) 



(151) 



When expressed in terms of the symmetric variables, the operator H 3 has been shown in 
Section ^| to leave invariant the subspace 

Mi = span{l,Ti,T 2 ,T 3 } . 

The restriction of if 3 to this subspace is explicitly given by 

3 3 

H 3 \mi= ~ Y} 2B + 1 + a ( 5 " M ^ T J ^ ~ C + ~ V + 1 ) T i+il ^ 

3=1 3=1 
3 3 

+ C + £}(4 - j) Tj .i 8 Tj - a J^(5 - j)(4 - 3)r 3 -2 d Tj + C +n + e . 

3=1 3=1 

The spectral problem for H 3 in the subspace M.\ thus reduces to diagonalizing the 4x4 
matrix of H 3 \ Ml with respect to the basis {1, r 1? r 2 , r 3 }, given by 



/e 


3C+ 


—6a 


\ 


^+ 


e — /5 — 2a 


2C+ 


-2a 





2C+ 


e - 2/5 - 2a 


c+ 


\o 





3C+ 


eo - 3/3/ 



with 

/5 = 2(a + 5) + 1. 

The eigenf unctions of i7 3 |yn 1 are 

Xj ( r ) = k j0 + kji t x + k j2 T 2 + T 3 , < j < 3 , 

where 

1 



[(£j - 2a + 2/5) (^ + 3/5) - 3C 2 ] 



**i = 7^2 + 2a + 2/5) + 3/5) - 3C 2 ] 



Ji 6C 2 
_ g J+ 3/5 

^ 2 " IcT 
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and Sj is one of the four roots of the characteristic polynomial 

£ 4 + 2(2a + 3/?)£ 3 + (4a 2 + 18a/? + ll/? 2 - 10C 2 )£ 2 

+ 6(/? 3 + 3a/? 2 + 2a 2 /? + 2aC 2 - 5/?C 2 )£ + 9C 2 (-2/? 2 + 2a/? + C\ ) . 

The corresponding eigenfunctions ^(x) = yu(z)xj(r) (0 < j < 3) of if 3 have energy 
E 



j — Sj + €q. If the parameters a, i? and C+ take the values 

2 , B = — , C_i_ = — 1 , 



a 



2 



compatible with the normalizability of the algebraic eigenfunctions, then the eigenvalues 
can be calculated numerically: 

E = -69.7926 
Ex = -64.1121 
E 2 = -56.4954 
E 3 = -44.5999 . 

The gauge part of the algebraic eigenfunctions is given by 

Xo{t) = 0.211595 + 0.376202 t x + 0.347522 r 2 + r 3 
Xl ( r ) = -0.959207- 0.00694338 n - 1.54596 r 2 + r 3 
X2 (r) = 0.00509379 + 16.3594 n - 4.08486 r 2 + r 3 
X3 ( T ) = -967.257 + 80.6047 n - 8.05004 r 2 + r 3 . 

Since the symmetric variables 

Ti = cosh xi + cosh x 2 + cosh x 3 

r 2 = cosh Xi cosh x 2 + cosh Xi cosh x 3 + cosh x 2 cosh x 3 
r 3 = coshxi coshx 2 coshx 3 

are all positive, xo, and therefore ?/>o, does not vanish in the region (pf). This shows that 
^0 is the ground state of the potential (|149|) . 

In the QES cases we are able to compute a few energies and their corresponding 
eigenfunctions algebraically, but we know nothing about the rest of the spectrum. In 
fact, it is natural to expect that the potential of this example has some levels between 
the ones we have obtained. 



5.3 ES potential from Case 3 

For the ES potential (B3|) 



E 



-2xi 



92 

4 



csch^ 



(152) 
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the gauge Hamiltonian in the symmetric variables is given by 



H 



N 



i.j L k=l 

Y^j[2B + l + a(2N-j-l)]r j d T . 
j 

C- y %2(N-j + l)T j - 1 d T .+E , 



d T d T . 



(153) 



with 



B 



(154) 



and E given by ( |131| ). The gauge Hamiltonian ( |153| ) preserves each subspace Ai k for 
arbitrary k. For instance, the restriction of Hn to M.2 has the following matrix with 
respect to the basis {1, n, T2, r\ } of .M2: 



/£ -iVC_ \ 

£o + l-/3 -(iV-l)CL -2iVC_ 

E + 2(a-(3 + 1) 4 

y E -2pJ 



(155) 



where /3 is defined by Eqs. (|141|) and (|1 54|) . By diagonalizing this matrix, we can again 
exactly compute four eigenfunctions with their corresponding energies for the multipar- 
ticle potential ( |1 52j ) . The energies can be simply read off the diagonal elements of the 
matrix ( |155| ), and coincide with those of the previous example (with B now given by 
( |154| )). The corresponding eigenfunctions are given by 



< % < 3 



where 



/1 = Y\ si nri ' 

j<k 



a / 3 



Xj X k 



)' CXP {? 



^(C + l + a(N-l))x k -^e- Xk 



and 



Xo 



Xi = n + 



2 2r 2 2CJNa + l) N(aN + l)C 2 

X2 = r x — + - — , tttt; , ^ r i + 



a + 1 (a + 1)09 + 1) 



P(a + l)(J3 + l) 



CJ1-N) 

X3 = T 2 + ~ — T T l 



N(N - l)Cl 



2a -(3 + 1 2(a-/3+ l)(2a - (3 + 1) 



(156) 
(157) 

(158) 
(159) 

(160) 

(161) 
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with 



k 



r-2 



j<k 



(162) 



Again, ip — fi has no zeros in the configuration space ([470 , an d is therefore the ground 
state of the system. As in the previous example, we have exactly computed some eigen- 
values and eigenfunctions of the potential ( |152|) for an arbitrary number of particles. 

5.4 QES potential from Case 3 

As a final example, consider the QES potential (|90|) with m = 1 



V N (x)=J2 



2 / Xj Xfe 



(163) 



The restriction of the gauge Hamiltonian to the subspace 

Mi = span{l,T!, ...,r N } 
is given in terms of the symmetric variables Tj by the following expression 

H n \mi = eo + C + n f 1 _ S r J ^ + ^ + ^ r i+ 1 ^ 

- ^ [C + a(2iV - j - 1) + 2] j Tj d Tj - C_ J2( N -3 + 1) tj-! d Tj , (164) 



with 



N 



(C + a (N - 1) + l) 2 + \ (N 2 - 1) + 2 



(165) 



The matrix (hij) 0<i - <N of if A r|^ 1 with respect to the basis {1 = r , r 1; . . . , rjy} of Mi is 
therefore tridiagonal, with nontrivial elements 



h jj-i = jC+ ) fyj+i = -C- (-W - j) , 
fyj = eo - j [C + a(2iV — j — 1) + 2] = . 



(166) 
(167) 



The tridiagonal character of the matrix (/i^) has important consequences for the calcu- 
lation of the eigenfunctions of lying in .Mi. Indeed, if we write one such eigenfunction 
with eigenvalue E as 



N 
3=0 



T 



3 1 



168) 
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it follows from Eqs. (|166|) -( |i~67j ) that the coefficients itj{E) satisfy the following relations 

jC + 7Tj-i(E) + {h 3 - E)^(E) -C.(N- j)n j+1 (E) = , j = 0, 1, 2, . . . , N , (169) 

where we can take 7r_i = and, without loss of generality, tt = 1. Let us now regard 
the energy E in Eq. (|169|) as a real variable, and let us consider the recurrence relation 



( |169| ) for arbitrary values of j G N. Following the usual procedure, we introduce new 
normalized coefficients ^j(E) through the formula 



*iW = Tlvr ~7\7 (-C-yn 3 (E) . (170) 



Note that ttq(E) = 1 implies that ttq(E) = 1. The coefficients iTj{E) then satisfy the 
following three-term recurrence relation in canonical form 

7T j+1 (E) — (E — hj) TTj (E) - a 3 iTj-i(E) , (171) 

with 

aj = -j(N-j + l)C-.C+. (172) 

The functions 7fj{E) defined by the latter relation with the initial conditions tc^i(E) = 0, 
ttq(E) = 1 are monic polynomials of degree j. It is well known that the canonical form of 
the recurrence relation ( |172j ) entails that the polynomials TCj(E) are an orthogonal family 
with respect to a suitably defined weight functional, [0. Furthermore, the fact that the 



coefficient aj vanishes for j = N + 1 implies that this functional is not positive definite; in 
fact, the polynomials iij with j > N + 1 must have zero norm. Therefore, the polynomial 
family {n j(E)} . >Q is weakly orthogonal, cf. WA - We have thus associated in a natural way 



a weakly orthogonal polynomial family to the QES many-body potential ( |163| ). Recall, 
in this respect, that it is possible to construct a weakly orthogonal polynomial family for 
(almost) every QES one-particle potential on the line, [18[ . 



Going back to the calculation of the eigenf unctions of H n belonging to .Mi, note first 
of all that their energies are the zeros of the critical polynomial 7Tjv+i (-£?). This can be 
seen, for instance, by observing that the j-th principal minor 5j(E) of the matrix E — (hij) 
satisfies the same recurrence relation as 7Tj(E), with the same initial conditions, so that 



and, in particular, 



8 j {E)=ix {E), (173) 



5 N+1 (E) = det(E - (hij)) = tt n+1 (E) . (174) 



Since C + C_ < for the potential ( |163| ), it follows from (|172|) that ctj is positive for 
j = 1,2, ...,N. This implies, [0, that the critical polynomial itn+i(E) has N + 1 
distinct real roots Eq < E± < • • • < En- By the previous equation, the spectrum of the 
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restriction of H N to M.\ consists of the N + 1 real eigenvalues E < ■ • • < E N . The 
respective (unnormalized) eigenfunctions x% = XEi are given by 

1 _^ . (N — j)\ 

where the coefficients TTj{Ei) are calculated either from the recurrence relation ( |171|) (for 
= or by computing the minors Sj(Ei). The corresponding eigenfunctions of the 
physical Hamiltonian are obtained by multiplying each function \% in ( |175| ) by the 
gauge factor 



j<k 



J] sinh " . exp J 1 £ [(c- + a(iV _ 1) + l) Xfc + C7 + eF* - C.e^] 1 . 

i<fe ^ ' I fc J 



(176) 



We can also express the eigenvalues of H N \ Ml and their eigenfunctions using continued 
fractions, by working with the quotients 



Then Eq. ( |171|) becomes 



q,(E) = -^§-, J = 1,2,.... (177) 



ry ■ 

q j+1 (E) = E - hj - -f- (178) 
Qj{E) 



from which we immediately obtain the continued fraction expansion 

q J+1 (E) - E - h, - E _ E _ hj _ 2 _ ■■■Y^h^E—h, (179) 

The eigenvalue equation in this formalism can be obtained by imposing the vanishing of 
qN+i(E) (since it is proportional to ttn+i(E), by Eq. (|177|) ), that is 

E-h N = "f- 1 ^ (180) 

E — tiN—i — E — hj\[—2 — E — hi — E — ho 



Alternatively, solving Eq. (|178|) for qj we obtain 



and an equivalent form of the eigenvalue equation ( |180|) follows by expressing qi(E) = 
E — h as a continued fraction: 

E — h 1 — E-ha- E-h N . x - E~h^ (182) 
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(since q^+x^E) vanishes if and only if E is an eigenvalue). The eigenf unctions are still 
given by ( |175|) , where of course now 

j 

n J (E) = Hq l (E). (183) 
i=i 

Consider, as an example, the potential ( |163| ) with N = 3 particles, for which 

e = -- A (C a + 2a + lf -2a 2 ~\C-C+. (184) 

The eigenfunctions of H 3 \m 1 can be expressed as 

Xj{r) = 1 + kj! n + k j2 t 2 + k j3 t 3 , (185) 
where the coefficients kji = TTi(Ej) are given by 

kfl = - g§r . (186) 

k * = ^2 i S J + ( 2a + P) e i + 3C - C +] > ( 187 ) 

k j3 = [£j + (4a + 3/3)f ? + (4a 2 + Qa/3 + 2(3 2 + 7C_C+)£ j + 6(a + /3)C_C+] . 

(188) 

In the latter equations we have set 

P = 2a + C + 2, (189) 
and Sj = Ej — eo is one of the roots of the polynomial 

tt 4 (£ + e „) =S 4 + (4a + 6(3)S 3 + (4a 2 + 18a/? + 11/3 2 + 10C_C+) £ 2 

+ 6 (f3(2a 2 + 3a(3 + (3 2 ) + (2a + 5/3)C_C+) S + 9C_C+ (2/?(a + /?) + . (190) 

For instance, if the parameters in the potential are given by 

a = 2, C_ = 2, C Q = 1, C+ = -l (191) 



then 
and 



eo = -32 

E = -54.5584 
Ex = -49.5497 
E 2 = -42.4323 
E 3 = -31.4596. 
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The eigenf unctions of if 3 in Aii are 

Xo ( T ) = 1 + 3.75974 n + 10.6142r 2 + 20.4323 r 3 
Xi(r) = 1 + 2.92495 ri + 4.5394 r 2 -3.94697 r 3 
X2 (r) = 1 + 1.73871 n - 0.496778 r 2 + 0.141027 r 3 
X3 (r) = 1 - 0.090071 Ti + 0.00986437 r 2 - 0.00137384 r 3 . 

Since the symmetric variables, given in this case by ( |162| ) and 

r 3 = z Xi+X]+Xk , (192) 

i<j<k 

are all positive, the function i[) = uxo is once again the ground state of the Hamiltonian 
H 3 . 



6 Summary and Conclusions 

We have used Calogero's construction of classical solvable many-body systems and applied 
it to the most general one-dimensional quasi-exactly solvable normalizable Schrodinger 
operator in the line. The corresponding quantum many-body Hamiltonians have been 
shown to have an algebraic structure that enables us to compute part or all of their 
spectrum by straightforward algebraic means. In all cases, if the one-dimensional seed 
potential is ES (QES) then the corresponding many-body potential constructed on it is 
also ES (QES). 

The QES examples include a generalization of the sextic deformation of Olshanetsky- 
Perelomov's rational model recently found in ]l(J, as well as some new deformations 



of the standard hyperbolic BCn potential. We have also found an additonal QES defor- 
mation of the .Btv rational model not related to the one considered in Ref. ]10| . 

The ES cases include Ajq and rational models with harmonic force and some 
hyperbolic OP potentials of BCn type, as well as the A N hyperbolic model with an 
external Morse potential discussed by Inozemtsev and Meshcheryakov, [fUjl - Since we have 
only looked for Hamiltonians having normalizable eigenf unctions, the models treated in 
H with no discrete spectrum do not appear in this work. The fact that we have considered 
only non-periodic one-dimensional seed potentials, |[20|| , explains why the trigonometric 
integrable Hamiltonians such as Sutherland's model are also absent. In a future work we 
plan to perform a similar analysis using periodic seed potentials which could lead to new 
results on the elliptic case, where much less is known. 

The ES problems have the additional feature that their associated gauge Hamiltonian 
H N preserves not only the subspaces M.^ but the smaller subspaces M.^ defined by (|34|). 
This implies that the number of particles N represents no difficulty in these cases, and 
allows us to calculate some levels and eigenfunctions for arbitrary N. Moreover, we have 
shown that an ordered basis in each subspace M.k can be chosen so that the action of 
Hn is triangular. As a consequence, the diagonal terms automatically give the discrete 
spectrum, and recursive expressions for the eigenfunctions can in principle be written. 
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It should be stressed that the results we have obtained depend greatly on the ansatz 
(j37l) for the gauge factor /i. Our ansatz has been chosen by analogy with the one- 
dimensional case, but it is clear that more general ansatze could lead to more general 



results. In his recent work , Calogero investigates a very general ansatz which gives 
rise to higher-body interactions. A systematic analysis of this method should focus on 
the equivalence problem, attempting to solve it in full generality. 
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8 Appendix 

The following expressions are useful in the transition from the canonical coordinates z to 
the symmetric variables r. (All indices run from 1 to N, with the restrictions indicated 
under the summation symbol.) 



Ai W - E 



4 



Mi Zk Zj 



Af ] = (193) 
aS 1} = ^N(N - 1) (194) 
Ap } = {N- l)n (195) 



(p) _ \^ 4 



A-2 = E 



i#fc#w (**-*i)(**-*») 



= A« = o (196) 
A^ 2) = -N(N- l)(N-2) (197) 



The following expressions are useful to express differential operators in terms of the 
symmetric variables (it is understood that = for k > N + 1): 
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(0) 



(N-j + l)T. 



3-1 



5 



(i) 
i 

>(2) 



*f = nrj - (j + l)r j+1 



J'-l 

^S? = (iV - i + l)ri_irj_i - y^(i - j + 2k)T i+k - 1 T j - k -i 

k=i 

3-1 

A if = £(* -J+ 2k + l)T i+ kTj-k-l 
k=0 

3 

4? = 3 r iTj + £(j - i - 2k)T i+k Tj- k 

k=l 



C7j°) = -(JV-j + l)(iV-j + 2)r j _ 2 
C7f = (JV-j)(JV-J + 1)^-1 
Cf =j(2N-j-l) Tj 
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